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^ ' It has been shown that there is a small possibility to experimentally violate 

the CHSH Bell inequality in a "classical" context. The probability of such 
a violation has been estimated in the framework of a classical probabilistic 
model in the language of a random-walk representation. 
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G^ . The Bell inequality (BI) is supposed to discriminate between "classical" world and 

"quantum" one. Any classical theory, i.e. causal theory governed by local (possibly, 
1^ _ hidden) variables, should fulfill the BI. Furthermore, any possible violation of the BI is 
Q ■ usually interpreted as a sign of quantumness. In other words, it is commonly believed 
^ ■ that only quantum mechanics is allowed to violate the BI. Then, the important question 
emerges: is it possible to violate the BI in the framework of a classical theory? If so, the 
powerful role of the BI could be a little diminished. It appears that the answer can be 
^ \ in principle positive, but it depends on details. Actually, the first such case is related to 
■ the notion of the postselection [Ij , other one refers to the notion of the memory loophole 
[2], there is an analog of the detection loophole presented in the NMR context in [3], 
and finally in optics in [4j. One should emphasize that a lot of them can be observed in 
laboratories. 

In this letter, we will show, and this is our main aim, that in finite number of mea- 
surement rounds it is possible to "classically" violate the BI, and we will estimate the 
probability of such a violation in the framework of a simple classical probabilistic model. 
To make things as simple as possible, we will confine ourselves to the BI in the version 
of Clauser, Horne, Shimony and Holt (CHSH) (see, [5] for a contemporary introduction), 
and a simple classical Bernoulli-like model. 
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We would like to strongly stress that our considerations have nothing to do (though, 
perhaps, could be somehow linked to) with real errors inevitably being encountered in 
actual measurements of true experiments (see, our toy model as a simple demonstration 
of the point). Otherwise, our conclusions would be entirely trivial. In other words, our 
"measurements" are exact, but our model is statistical. 

First of all, we would like to remove the seeming contradiction between the well-known 
proof of the CHSH inequality and our a bit controversially sounding statements. Namely, 
according to our interpretation, which plays an auxiliary role and is not crucial for our 
model, one should note that a delicate point in the CHSH argumentation consists in 
summing up (with one minus sign) and collecting the all four terms ai{\)bj{\), i,j = 1, 2, 
i.e. 

ai(A)6i(A) - ai(A)62(A) + a2(A)6i(A) + a2(A)62(A), (1) 
under the common probability measure P(A)dA, yielding 

J P(A) [ai(A)6i(A) - ai(A)62(A) + a2(A)6i(A) + a2(A)62(A)] dA. (2) 

Here, as usual, aj(A) = ±1, bj{X) = ±1 are values of the polarizations measured by the two 
(distant) observers A and B, respectively. The indices i,j = 1, 2 label the two orthogonal 
directions of the polarizers A and B, and A is a collective hidden variable governed by the 
probability distribution P{X) (see, [6] for a very short introduction). 

According to our interpretation the loophole in the CHSH argumentation consists in 
the direct addition and integration in ([T]) and One should take into account that, 
simply speaking, each of the four correlations in is, in practice, being measured in 
different time instants. Since hidden variables are in principle allowed to evolve in time, 
we have in general different A's in each of the four terms, and ([T]) and ([2]) does not, in 
general, make sense. Seemingly, similar doubts, but formulated rather in mathematical 
than physical terms are presented in [U El IS] . 

If somebody is still sceptical concerning our above interpretation with the time t in its 
central role we offer the following additional remarks: 1) Let A = X{t) = t, for example. 
Does Eq.(([l])), Eq.((l2])) or the standard normalization J P(A)dA = 1 make sense? Rather 
not. 2) The role of time is to break correlations, which forbids making the valid sums 
of the type Oi ± 02 being used in the proof of the CHSH inequality. Our mathematical 
colleagues have already noticed some inconsistences of this type, referring it to the issue 
of different probabilistic spaces (see, earlier references). 3) Our toy model presents a kind 
of direct "experimental" explanation of this phenomenon. 4) Finally, one can skip this 
explanatory paragraph at all, if one dislikes it, and directly jump to our model and the 
results which are interpretation independent. Our interpretation plays only an auxiliary 
role in our paper. 

Obviously, the very fact that the proof has a loophole does not automatically mean 
that the CHSH inequality can be actually classically violated unless we show it explicitly. 
That will be done in due course. 

There is also another seeming contradiction in our proposal. Namely, nobody is re- 
porting an experimental classical violation of the CHSH inequality due to finite statistics. 
The explanation is very simple. Actually, the probability of such a violation in a real 
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experiment is very small. In fact, it depends on the number of the measurement rounds. 
It will be shown that the greater the number of the measurement rounds the lower prob- 
ability of the violation. The fact that the number of measurement rounds is finite is 
crucial. 

A conclusive explanation of the both above mentioned seeming contradictions will 
be given in the form of the proposed model. Mathematically, the model is formulated in 
probabilistic terms. In other words, formally, speaking in the language of hidden variables, 
one could state that the hidden variables are purely random entities. From physical point 
of view, we could interpret the randomness of the hidden variables in various ways. They 
could be fundamentally random. Their randomness could follow from complicated internal 
classical statistical mechanics. Or possibly, their classical mechanical evolution could be 
fundamentally very complicated, e.g. chaotic. Anyway, we are only interested in this letter 
in a principal possibility of classical violation of the CHSH inequality independently of a 
possible physical mechanism, i.e. the nature of the hidden variables. 

A notable precedent of our work is a multi-thread and polemical study undertaken 
by Gill in [TU] and [TT]. One of the relevant threads of the Gill's papers is a polemic 
centered around the role of statistics in violation of the CHSH Bell inequality. Namely, 
some authors propose classical scenarios significantly, according to them, violating the 
CHSH Bell inequality. Gill argues, using advanced statistical calculus, that such large 
violations, suggested by them, are impossible. 

To quickly illustrate the point, first, we will present a toy version of our model. To 
this end, let us confine ourselves (in this toy version) to only four measurement rounds 
(in the full version that number will be large). Let us perform the four measurements in 
four different time instants, denoted: 1, 2, 3 and 4. We could even identify time with the 
hidden variable governing the process, but it is optional. "Accidentally" , it appears that 
our experimental data are such as those given in Table 1. Obviously, nothing forbids to 
obtain such data. As it is clear from Table 1 we deal with maximal possible violation of 
the CHSH inequality. Namely, we obtain the (maximal) number 4 instead of the classical 
bound given by the number 2 or the quantum bound given by 2\/2. We can observe, 
enriching our earlier arguments, that the constraints used in the standard proof of the 
CHSH inequality do not work because the variables a and b in different time instants 
(or for different hidden variables) are totally unconstrained. Evidently, the data given in 
Table 1 constitute 
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Table 1 : Data maximally violating the CHSH BI in a classically allowed experiment 

a very particular configuration out of many others. How many? A simple calculus, 
reduced for simplicity to c = ab (which is equivalent to a more tedious calculus given in 
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terms of a and b separately) gives 2^ = 16 all configurations. We have assumed that the 
i,j pairs are fixed. This is a simplification (also present in our full model) — in a more 
realistic situation only the total number of the measurement rounds should be fixed (to 
4, in this example). Since there are 2 violating configurations (one for 4 and one for —4), 
the total probability p of the violation of the CHSH inequality in our toy model amounts 
to 

2 1 

p= — = - = 0.125. 3 

^ 24 8 ^ ^ 

Now, let us consider the full model. In other words, we will estimate the probability 
of the violation of the CHSH inequality, in the spirit of the toy model, for a larger number 
of measurement rounds. In terms of probability theory our model is described by a 
Bernoulli process [12]. But it is more convenient, from intuitive point of view as well as 
from calculational one, to recast the model into (four-dimensional) random walks. As 
we have just observed, we can think and work directly in terms of the outcome c = ab 
instead of in terms of a and b, separately. In our random-walk representation the value 
c = +1 corresponds to the step forward, and the value c = — 1 corresponds to the step 
backward. The direction forward/backward (±1) is "decided" first. Next, the "particle" 
(we mean the abstract particle of the random-walk representation) should be "informed" 
which dimension out of 4 should be followed. There are 4 dimensions corresponding to 
the 4 possibilities (pairs) given by the orientations of the polarizers A and B {i = 1,2 
and j = 1,2). The total number of the steps in one of the four dimensions is denoted by 
rii ii = li2,3,4). Since each step assumes the value -|-1 or —1, the actual ith coordinate 
of the position of the "particle" is the numerator of the zth correlation. In other words, 
each of the 4 correlations entering to the full CHSH correlation is of the form ( "frequency 
definition") 

^ = ^. (4) 

Hi Hi 

Therefore, the experimentally read CHSH correlation is given by 

ni n2 ns n4 ' 

Now, the inequality corresponding to the (rare) violating cases we are interested in as- 
sumes the form 

\C\ > 2. (6) 

For a large number of the measurement rounds, i.e. 

rii > 1, (7) 

we can use a continuous approximation. The probability distribution of the one- dimensional 
discrete random walk fL2\ is summarized by: 
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Fig. 1: Probability of the violation of the CHSH BI versus the total number 
of the measurement rounds N {N = ni+n2 + ns + n4), in a doubly-logarithmic 
scale. 3 variants — solid line: rii — n2 = = n^; dashed line: lOni = 712 = 
713 — ^4/ dotted line: lOOni — 712 — ns — n^. The five vertical intervals 
correspond to the exact (discrete) calculations easily performed for small N 
and equal Ui. The lowest points of the intervals correspond to C > 2, whereas 
the highest ones correspond to C ^ 2. 

each coordinate is independent. Therefore, the full four-dimensional probability measure 
is given by 

1 X 1 1 ^l'^^2^ 3:4^ 

dP„i„2„3„4(a:i,X2,X3,X4) = -j==-^==-^==-== e'^'^'^-s'^ d^x. (10) 

^/ZTTni y/2TTn2 \J 2'Kn-j, \J 21^71^ 
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From mathematical point of view the task is to calculate the probability of finding the 
walking particle in the four-dimensional space outside the layer L bounded by the two 
hyperplanes C2 and C_2 (see, Eq. ([5]) and ([6])), where 

C±2 : - - - + - + - = ±2. (11) 

ni 

For technical reasons we will "isotropize" the coordinates of the probability measure by 
the following change of variables: 

Xi = ^/2^iZ,, z= 1,2,3,4. (12) 

Now, the ( "isotropic" ) probability measure ffTOj) is of the form 



'^Pnin2nzni{.^l, Z2, Z^, Zi) — — e ^* d'^Z, (13) 



whereas the equations of the bounding hyperplanes (ITTI) are 

C±2 : J— zi - J— Z2 + \ — z^ + \ — Zi = ±2. (14) 
\ ni \ n2 V ^3 V "-4 

The distance d between the hyperplane given by the equation ^ ajX, = b and the begin- 

i 

ning of the coordinate system is expressed by the formula: 

'2^ at 



Making use of the isotropy of the probability measure we can rotate the hyperplanes so 
that they are orthogonal to the coordinate Zi, say. Therefore, the probability p we are 
interested in is given by the integral 



oo 



p(7li,7l2,n3,7l4) =2-^ j (\Zi j d' 

— oo 

dt e-*' = erfc(d) (16) 



d 

oo 



TT 

d 



=erfc 



- + ^ + ^ + ^ ' ' 

ni 712 ria "4 



where in the last line we have made use of 0141) and 0151) . Since p is very quickly damped 
by large Tij, the probability of the violation of the CHSH inequality for ?ij 1 is very 
small. 

A visual summary of our results is given in Fig. 1. Violation of the CHSH BI, as well 
as its scale, is now evident. For simplicity, in our model, we have assumed fixed values 
of Tij. In a more realistic model only the sum N = ni + n2 + + n^^ should be fixed. 
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which would roughly correspond to a mean of ( fT6ll with respect to rij. But, obviously, 
our technical simplification, qualitatively, does not change our final conclusion. 

The continuous approximation is only a simplifying, technical trick because discrete 
calculations, as combinatorial ones, are straightforward only for small A^. Nevertheless, 
we observe that the continuous approach does work quite good also for small A^. 

Not to cause any misunderstandings in the field we would like to stress that in the 
limit of infinite number of measurement rounds the CHSH Bell inequality is restored in 
the classical domain. In other words, assuming the "infinity" condition as a premise in the 
CHSH Bell theorem removes the classical violation but evidently the "infinity" condition 
is experimentally unrealistic. 
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